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Abstract 



We investigate the equilibrium properties of arrays of coherent strained 
islands in heteroepitaxial thin films of bidimensional materials. The model 
we use takes into account only three essential ingredients : surface energies, 
elastic energies of the film and of the substrate and interaction energies be- 
tween islands via the substrate. Using numerical simulations for a simple 
Lennard-Jones solid, we can assess the validity of the analytical expressions 
used to describe each of these contributions. A simple analytical expression 
is obtained for the total energy of the system. Minimizing this energy, we 
show that arrays of coherent islands can exist as stable configurations. Even 
in this simple approach, the quantitative results turn out to be very sensitive 
to some details of the surface energy. 

Keywords : Semi-empirical models and model calculations, Epitaxy, Self-assembly. 

I. INTRODUCTION 

Manufacturing arrays of three dimensional nanoscale islands is a field of intense research, 
with many potential optoelectronic applications ('quantum dots') The production of 
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these structure represents a difficult challenge. One possible way to proceed is by spon- 
taneous formation of coherent (i.e. dislocation free) strained islands in heteroepitaxy : 
this approach represents an elegant and simple way to produce quantum dots. For instance, 
when growing a layer of InGaAs on GaAs(OOl) or of SiGe on Si(OOl) spontaneous 
formation of three dimensional coherent islands is observed after depositing a few layers. 
However, our understanding of the mechanisms that lead to island formation from strained 
layers is still incomplete, and even the equilibrium properties of such structures are not well 
understood 

A fiat strained layer is unstable P-|lO|] against island formation. The driving force of this 
instability is the elastic energy gain due to the relaxation of the stress at the top of the 
islands. Basically, the elastic energy gain for the creation of islands is proportional to the 
volume V of the islands (per unit area of an island). The surface energy loss, on the other 
hand, scales with the volume as V"^^^. Hence it is clear that, above some critical volume, the 
growth of the island will be energetically favored. In this simple description, the equilibrium 
state of the system consists in a single very big island. 

Experiments, however, reveal a different behavior. Regularly spaced islands with a nar- 



row size distribution |TT| are observed on the substrate, sometimes above a thin wetting 
layer. To explain the difference between the simple argument and the observations, several 
authors have suggested that the density of islands is fixed by the kinetics of the growth pro- 



cess 



T^JT3|| . However, some systems of islands are not affected by reasonable annealing [jT4[ 



Either the driving force becomes very small so that the system never reaches its equilibrium 
state, or there is a stable or a metastable state corresponding to a given island density. 
Indeed, Shchukin et al. [|1^ have argued that coherent arrays of tridimensional islands may 



be thermodynamically stable, when taking into account elastic energy, facets energy, surface 
stress and the contribution of the edges. 

In this paper, we study the equilibrium properties of a bidimensional strained layer. 
Our aim is to gain an understanding of the structure by using a minimal model, with 
the simplest possible ingredients. To describe the energetics of the layer, we only take into 



account surface energies, interactions between islands, and the elastic energy of the substrate 
and the strained layer. We have also considered the possibility that the islands grow on a 
wetted surface so that our analysis is also relevant to the Stranski-Krastanov growth mode. 
Some of the terms that contribute to the energy cannot be computed analytically, or only 
with uncontrolled approximations. When this is the case, we use simulation of a simple 
model (Lennard- Jones interactions) to investigate the approximations done in the analytical 
calculation. 

Minimizing the total energy of an array of islands, we find that, using only these simple 
ingredients, stable arrays of islands are predicted in a range of parameter values. The size 
and density of the islands at equilibrium and the height of the wetting layer are also obtained. 

In the first part of the paper, we define the system we study, and the model used in 
the simulations. In the second part, we will precisely review all the terms that enter the 
energy calculation. Elastic terms are calculated using linear elasticity theory, supplemented 
by atomic simulation results. Results are presented in section |rV|, and briefly discussed in 
section |V] 



II. DEFINITION OF THE SYSTEM AND SIMULATIONS. 

A. The system 

In this paper we compare the energies of a strained layer containing an array of coherent 
islands, with that a fiat strained layer of the same volume. For the sake of simplicity, we 
restrict ourselves to a bidimensional geometry. In figure |l|, the two situations for which 
we want to compare the energies are described schematically, and the geometric parameters 
relevant to describe each configuration are introduced. In both situations (uniformly strained 
layer or regular array of islands), 9 is the amount of deposited particles per unit of length. 
9 is also the height of the uniform layer in figure 0(a). In figure 0(b), h and i are the height 
and width of the islands P,|16[ so that the volume of an island is = M. We choose the 
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angle $ equal to 7r/3. z is the height of the wetting layer. The distance between the islands 
is L^. Note that an infinite value of corresponds to a substrate with a single island. Since 
the amount of matter deposited is the same in both configurations (a) and (b), we have 

ei^ = zL^ + hi (1) 

Our aim is to calculate the difference of energy AE between the two configurations shown 
in Fig. H. AE is a function of i, h, Lx, z and 6. 6, which is the amount of adatoms deposited 
on the surface, is a controlled, fixed parameter. Since Eq. |l] connects the four variables 
i, h, Lx and z^ AE is a function of only three independent variables. We choose V, r and z 
as independent variables where, V = M is the volume of an island and r = hjl its aspect 
ratio. 

By minimizing the difference of energy per particle AE/OLj. at constant 9 and for a given 
V with respect to z and r, we obtain AE /OL^ as a function of the volume V of an island at 
given coverage. We can then locate the equilibrium states by minimizing with respect to V . 



B. Simulations 



To evaluate some elastic terms in the difference of energy AE or to check some as- 
sumptions in our analytical calculations, we have performed simulations with a 'Lennard 
Jones material'. In this material, particles interact via a Lennard- Jones potential. The 
Lennard- Jones interaction potential is given by: 



V{r) = 4e 



(T\i2 fa 



(2) 



Again for the sake of simplicity, we consider only first neighbors interactions. This implies 
that the strain at the surface is identical to the strain present in the bulk material (no sur- 
face relaxation), so that there is no excess strain or stress associated to the surface. These 
model materials therefore have no surface stress. We use a bidimensional description and 
assume that the strain free material (substrate and layer) has a triangular lattice structure. 



Since there are two types of particles, 3 couples of parameters are needed to describe inter- 
actions between atoms in our system : substrate-substrate (e^s, Css), layer-layer (e^;, an) and 
substrate-layer (e^^, cr^s). We take the lattice spacing of the substrate as the unit of length, 
so that we have only to specify the misfit mo between the two materials. We assume that 
the value of aig is the average of ass and an. So that we haveQ: 

<^ss = ^ (3) 
du = assil - mo) (4) 

(^ss + O'li . , 

(yis = — ^ — (5) 

In our simulations, we choose the values of the bonding energies and of the misfit to be: 
e,, = 0.058e1/, = OMeV, en = OMeV and mo = 0.06. 

In order to limit the system size, we assume that only the atoms of the substrate in the 
first Usub layers close to the surface can relax. Since the Green function |T^ of a semi infinite 



system decreases as 1/r where r is the distance between an applied force on the surface and 
a given point in the system, the effect of limiting the depth of the substrate is basically to 
screen the horizontal forces beyond a distance Usub- Hence, two islands at a distance more 
than Usub do not interact in our simulations. We use a substrate of typical depth 60 (in 
substrate lattice spacing units) or more. We will check in part |III B| that at this distance, 
interaction forces between islands become indeed negligible. 

For a given topology (islands or strained layer), our simulation consists in relaxing all 
atomic positions towards the nearest energy minimum, using a standard conjugate gradient 
algorithm. We are then able to evaluate the elastic energy of the system by subtracting the 
sum of the equilibrium bonding energy for each pair of particles from the total energy of the 
system. 

Remembering that interactions are only present between nearest neighbors, we can com- 
pute the elastic constants and surface energies of the material (at zero temperature). The 
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surface energies of the substrate per surface atoms is 7^ = ess, that of the epitaxial layer 
7; = eii. The interface energy is 7^/ = (e^^ + — 26^; and the Young modulus of a bidimen- 
sional Lennard Jones crystal with a triangular lattice for the substrate Es = 
for the layer Ei = J^f^ ^^. We can also compute the Poisson ratio which is the same for 
both materials : u = l/2\^. With the numerical values we use as parameters, the strained 
material totally wets the substrate: 7^ — (7; + 7^;) > 0. 



III. INGREDIENTS 

In this section, we list the different contributions to the difference of energy AE be- 
tween the two configurations mentioned in Fig. |l[ Three essential ingredients are taken into 
account: elastic energies, interactions between islands and surface energies. Each of these 
contributions is discussed below. 



A. Relaxation energy 

The elastic energy per unit length of a film of height 6 with a uniform strain mo is given 
by the expression : 

El!:luJL. = ^mle (6) 

Where Ei is the Young modulus of the strained material. 

The exact expression of the elastic energy of a (substrate+island) system is, to our 
knowledge, unknown. Two studies ||T8|-pO| propose an approximate expression for this elastic 
energy, in the absence of a wetting layer {z = 0). In both studies, the elastic energy of an 
island is written as : 

= ^rnlR{r)V. (7) 

Here V is the volume of the island and R{r) is a dimensionless factor determined by the 
shape of the island and the elastic constants of both the substrate and layer materials. These 
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studies show that the shape factor depends essentially of the aspect ratio r of the island, and 
only weakly of the ratios between elastic constants. This shape factor takes into account 
both island and substrate contributions to the elastic energy. 

Fig. 1^ presents the shape factor R{r) calculated from our simulations for the peculiar 
case of islands with a vr/S contact angle. R{r) goes to 1 when r — : a very fiat island can 
not relax its elastic energy. As r increases, R{r) decreases, since atoms in the islands can 
relax as soon they are far enough from the substrate. In our simulation, we cannot evaluate 
the shape factor R{r) for aspect ratio greater than Vmax = tan(0) ^ 1.7 due to the choice of 
the islands shape. 



For the numerical calculations of Sect. IV, we need a accurate analytical expression of 



R{r). R{r) physically corresponds to the elastic energy relaxation due to the creation of 
an island, it will then appear in the only negative term of the expression of AE : our final 
result will thus be sensitive to this function. To obtain an analytical expression for R{r), 
we have fitted the results by an exponential : we find R{r) = 0.13 + 0.87 * exp(— r/0.18). 
We have checked that neither the expression of Ratsch and Zangwill (which contains a 
free parameter) nor the expression of Kern and Miiller P^JTP[ ] (no free parameter) allow to 



fit our simulation results (Fig. However, this disagreement may be due to the difference 
of contact angle between their islands and ours : Ratsch and Zangwill and Kern and 
MuUer islands have a contact angle tt /2 with the substrate, whereas our islands have a 7r/3 
one. Moreover, our simulation also show that the shape factor does not significantly change 
with the height of the wetting layer. 

B. Interaction between Islands 

To evaluate the interaction energies between islands, we use the same formalism, the 
same approximations and the same islands shapes (cf. Fig. [1|) as TersofT [|6|,p!6|. We assume 



that an island with a shape given by h{x), exerts on the substrate surface an horizontal force 
density given by = aod^h^x) where ctq is the bulk stress in a layer uniformly strained to fit 
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the substrate. We calculate the displacement of an atom of the substrate surface caused 



by the presence of an island at a distance using the Green function |jT^ of a semi infinite 



plane. We then evaluate the strain Uxx caused by the presence of the island : 



Ux = — r:^motan(</)) 



s 



{Lx + al2f - {l/2f 
{Lx - a/2)2 - (£/2)2 
{Lx + i/2y-{a/2Y 
{Lx - £/2)2 - (a/2)2 



+aln 



Uxx - g 



■'X 



where a is defined by a = h/ tan(0) and we have used do = Eirrio. i, h and are defined on 
Fig. m. Eg and Ei are respectively the Young modulus of the substrate and of the strained 
material. 

An island located at a distance Lx from another sees a misfit m = mo — Uxx different from the 
original misfit mo. From the elastic energy Eq. |^ of a single island, we deduce the interaction 
energy between two identical islands separated by a distance Lx'. 



Einteract — EiR{r)m'^ 



V (10) 



mo 

where Uxx is given by Eqs. ^ and ^ A first order asymptotic expansion of Eq. ^ shows 
that Uxx scales as —l/L^ for large Lx- Eq. |lO| then corresponds to a repulsive interaction 
between islands. Figure 01 shows a comparison between the interaction energy of two islands 



obtained from simulations using the Lennard Jones model and the result of equation |T0. 
Remembering that there is no adjustable parameter in Eq. 0, this expression leads to a 
reasonable estimate for the interaction energy. We have neglected the influence of the height 
of the wetting layer. A wetting layer, however can be taken into account in the simulation, 
and fig. |^ shows that it does not modify the order of magnitude of interaction energies. 

Moreover, the main effect of this interaction energy in the whole expression of AE 
consists in a strong short range repulsion. Thus, with or without the presence of a wetting 
layer, the approximate expression |TD| should allow to get the physical effects of the repulsive 
interactions between islands. 
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C. Surface Energy 



We have to calculate the surface energy difference between the two configurations of 
Fig. |I]. For a material without surface stress, and with the elastic terms calculated Sect, [ill A 



already taking into account the strain of the surface, we only have to consider the strain 



free surface (i.e. atoms are in their equilibrium position). Following Miiller and Kern ||22 



and simulations results of Wang et al. |T3|, the surface energy of the layer and the interface 



energy depend on the height of the layer, and we define S{z) : 

S{z)=^i{z)+^,i{z) (11) 
S{z) is the surface energy cost to create a wetting layer of height z on the substrate. When 



z tends to infinity, S{z oo) tends to 7^ +7^^ where 7^ and 7^; has been given in Sect. [II B 



When z tends to 0, the wetting layer disappears, and only the substrate is left, so that 



S{z — > 0) tends to 7^. As suggested by Miiller and Kern [|22|, we choose an exponential 
function for S{z) (valid for semiconductors and metals). This choice should qualitatively 
give the correct behavior for the function S{z). 

S{z) = (7, - 7i - 7,;) e-^/^° + 7i + 7s; (12) 

And following Miiller and Thomas |2^, we define surface energies as : 

7K^) = 7; (1 - e-^/^") (13) 
7sK^) = - 7.) (1 - e-^/^") + 7. (14) 

Where zq is a parameter comprised basically between 1 and 5. We have chosen 2:0 = 3 in 
this work. This choice of the screening distance zq, and the choice made for the function 
S{z) will be discussed in Sect. ^ We are now able to write the difference of surface energy 
between the two systems of Fig. |I|. 

^E^^rface = [S{z) - S{9)] [L, - i - a] 

+-fiiz + h) i-a + 2—^ 

sm[(p) 



+jsiiz + h) [i + a] 

-S{e)[e + a\ (15) 

With a = h/ tan(0). 

The first term on the right side of Eq. |15| is the surface energy difference for a film of length 
— i — a when changing the thickness from 6 to z. The last three terms correspond to the 
surface energy difference between a flat film (width i + a) and the island (height h and width 
i). We take = 7r/3 so that for a triangular lattice, island sides surfaces are equivalent to 
top surfaces. 



IV. RESULTS 



Using Eq. |10] and [T5|, we have all the ingredients to write the expression of AE. 

V 



AE = —ml 
2 ° 



+ [S{z) - S[Q)\ [Lx-i- a] 
+-fi{z + h) £-a + 2 ^ 
+-fsi{z + h)[i + a] 
-S{e)[i + a] 



sin{(f)) 



(16) 



Where Uxx is given by Eqs. | and ||. Using variables V, r and z, and writing the energy per 
particle, we obtain : 



AE E 



9L. 



-rrin 



i?(r)(l 



u 



XX \ 2 



- 1 



s{z) - s{e) 
e 



+-fi{z + VW) 



1 - ■ 

e-z 
ve 



e-z 
e 



V/r + Vrl^/tan(0) 

V 

\/ rV \/ rV 



tan(0) sin{(f)) 



+-fsiiz + VrV) 

e-z' 



V/r + 



'rV 



-siey 



ve 



V/r + 



tan(0) 

7v 



e-z 
ev 



tan(0) 



(17) 
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Where we have used h = VrV, £ = \jv/r and = from Eq. |1|. Each term of 
equation |1^ has been derived from the simulations. But we now assume that this equation 
is vahd for any value of the parameters and is independent of the details of the simulations 
(for example, the value of Ugub is not meaningful in equation [T7|) . 

For 9 and V given, we minimize Eq. |l^ compared to the variables r and z using conditions 
< z < 6 and {i + a) < (islands cannot overlap). 

Figure | presents the results of this minimization procedure: each curve corresponds to a 
fixed value of 9. For large enough coverages, a minimum in the curve is obtained at small 
volumes, corresponding to the stability of an array of islands. For greater volumes, the 
function ls.E/9Lx increases up to a local maximum, and then decreases up to a horizontal 
asymptote for infinite volume. This decrease is not shown on Fig. ^ because of the choice 
of the axis scale. In the case of Fig. ^, for each curve presented, the horizontal asymptote 
has a greater energy level than the local minimum at small volume, so that these minima 
correspond to a stable configuration. 

With the numerical values we use for the parameters, the flat layer is unstable for a 
coverage 9 larger than 1. The critical thickness is thus lower than 1. However, for 9 = 1, the 
volume of the stable islands is only 4 and the height and width is about 1.2 and 3.4. Thus, 
if these values were valid in an experience, the measured critical thickness would be greater 
than 1 since these islands would be of about the same size as the thermal perturbations of 
the fiat surface. For a coverage greater than 1, the energy is minimum for an island volume 
ranging from 4 to 1200 particles, depending on the coverage. Fig. ^ shows the evolution of 
the volume at the minimum of energy /S.E/9Lx as a function of the coverage 9. This volume 
is an increasing function of the coverage. At the minimum, we find that z = 0, i.e. there is 
no wetting layer. On Fig. ^, we also show the distance between islands in the minimum 
energy configuration. The island density decreases with the coverage 9. Fig. ^ shows the 
island height, width and aspect ratio as a function of the coverage at the minimum of energy. 
The aspect ratio remains small (lower than 0.7) for all the coverages. 
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V. DISCUSSION 



Our study shows that the existence of stable arrays of islands can be explained using 
relatively simple ingredients. The equilibrium state depends on the value of the coverage, in 
qualitative agreement with experimental observations ||^. Typical volumes for the islands 
are of the order of 500 particles, and typical widths are in the range from 20 to 40 (in 
substrate lattice spacing units). This corresponds to islands of width from 4.7 nm to 9.4 
nm, using the silicon lattice spacing. These widths are in good agreement with measured 
widths of "huts" in the Ge on Si (001) experiments [Q. 

Moreover, we find that there is no wetting layer at equilibrium. System such as Pd on 
Cu or InGaAs on InP |^ exhibit this type of behavior and even atoms of the substrate 



go in the islands. 

To investigate the sensitivity of our results to the choice of the function S{z) (equa- 



tion |T2D, we use the following definition of S{z) : 



S{z) = 7, - for 2 < 3 (18) 

Siz) =7/ + 7.z for^>3 (19) 

Fig. presents the result of the minimization of Eq. |1^ with this new expression of S{z). 
Fig. 1^ still shows minima but these are metastable and the energy barrier to leave them is 
smaller than in Fig. ^. We also find that our results are sensitive to the value of 2:0- 

Thus, comparing results of Fig. H and ^ could suggest that our model is very sensitive 
to the choice of the function S{z). Nevertheless, this sensitivity is mainly quantitative : the 
function AE/OL^ also shows a decrease at very large volume (not shown on Fig. H). There- 
fore, changing the function S{z) does not change the qualitative behavior of the function 
AE/6Lx : curves always present a local minimum at small volumes, a local maximum at 
intermediate values, and then a decrease to an horizontal asymptote. 

Concerning the quantitative disagreement between Fig. ^ and |^, our model requires the 
minimization of a function of several variables. The presence of smooth local minima of this 
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function may introduce a relative sensitivity to the input quantities. A plot of the functions 



described by equations 12 and 03 shows that the difference between these is not negligible. 



Finally, taking the derivative of expression O with respect to V, we note that the value of 



V minimizing Eq. depends on the value of dS{z)/dz. Unlike equation |T^, equation ^ 
has a zero derivative for z greater than 3 : this unphysical behavior may greatly affect the 
results of the minimization of AE/OLx- 

These results show that a very precise description of the surface energies would be needed 
in order to give the theoretical calculations a predictive character. This was not the aim 
of the present paper, but such a description may in principle, be obtained from ab-initio 



calculations 13 



Concerning the behavior of AE/OL^j. at large volumes, curves of Fig. |^ display an energy 
minimum for an infinite island volume for large enough coverages {6 > A). This minimum 
is the one predicted by the very simple argument of Sect. | : the corresponding state of the 
system consists in a single very big island. At this minimum, the height of the wetting layer 
is not but takes the value 3 for most coverages 6. As announced, the curves of Fig. ^ also 
show these minima at very large volumes, but the wetting layer has a height in the range 
from 1.3 to 1.7 depending on the coverage. Since the thermodynamical force that drives 
the system to this equilibrium configuration decreases to zero when the volume of islands 
V tends to infinity, it would not be surprising that the system does not evolve even after 
reasonable annealing because of kinetic limitations. Hence, our analysis suggests that some 
experimental observations of large islands are the result of kinetic factors and that some 
observed islands do not actually correspond to an energy minimum. 

Fig. ^ and |^ do not show any nucleation barrier for islands formation. This is actually a 
direct consequence of the minimization of AE / OL^ with respect to r and z : we are looking 
at the equilibrium states of a strained layer and a priori, these are different from the states 
occurring along the growth of the system. Basically, in our treatment, if the production of 
islands is energetically unfavorable, the equilibrium state of the system is the flat strained 
layer and the minimum of AE/OL^ is obtained for z = 9 {L^ has then an infinite value so 
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that equation |l] remains valid). Thus, by definition, reported values of Fig. ^ and |^ can not 
have positive values and no nucleation barrier can appear in these figures. Nevertheless, 
during the growth of a real strained layer, we expect that a nucleation barrier drives the 
islands formation. Indeed, before the formation of large islands, the system has to create 
smaller one which would not correspond to equilibrium states described by Fig. ^ and |^. 
The nucleation barrier depends on the kinetic pathway followed by the system to reach its 
equilibrium state : thus, the determination of the energy barrier demands a kinetic study of 
the transition. 

To check the robustness of our model, we performed few tests changing the values ei2 
and €22 of the bonding energies. We have changed the values of the ratio €12/611 from 0.9 
to 1.21 (this ratio was fixed at 1.03 previously) and the ratio e22/eii from 0.52 to 0.86 (0.66 
previously) : we find that the qualitative behavior is not affected. Curves always present 
minima at small volumes with no wetting layer, a local maxima and then a decrease up to 
an horizontal asymptote. 

The observations of Ref. shows that "huts" of Germanium on Silicon(OOl) disappear 
during the growth and simultaneously "domes" appear. Such result could be interpreted by 
an energy landscape such as the one of Fig. |^ : 'huts' would correspond to a metastable 
state whereas, 'domes' would correspond to a non equilibrium state. To definitely answer 
this question, one would have to describe very carefully the system Ge on Si(OOl). 

In Eq. |l^, we made the assumption that variables h, i and are continuous. This 
is obviously an approximation. The only alternative, however, would be a fully atomistic 
simulation. 

Our analysis does not take into account surface stress in the elasticity calculation in 
Sect. |T|. It is interesting that nontrivial stable configurations are nevertheless obtained, 
showing that surface stress is not an essential ingredient in island formation. Surface stress 
could be included automatically in the numerical simulations by taking into account long 
range interactions between atoms. Another theoretical challenge would be to take into 
account the possible formation of alloys during the growth : this again seems rather difficult 
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at the analytical level, since it involves taking into account new elastic coefficients, and new 
variables such as the volume of the alloy. 

Our work shows that, by combining an analytical approach with some ingredients ex- 
tracted from numerical simulations, a quantitative study of the stable configurations for a 
strained layer is possible. Although we have limited ourselves to a simple two dimensional 
geometry and a model Lennard- Jones potential, it seems possible to extend the approach to 
more realistic interactions (e.g. Tersoff ' potentials for Si and Ge) and a three dimen- 
sional geometry. 
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FIGURES 
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FIG. 1. a) Flat strained layer, b) Strained layer containing islands. 
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FIG. 2. Shape factor R{r) as a function of the aspect ratio r. The solid Hne has been obtained 
from the equation : R{r) = 0.13 + 0.87 * exp(— r/0.18). We have also compared our simulations 
with the predictions of Kern and Miiller and Ratsch and Zangwill. 
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FIG. 3. Interaction Energy between two islands as a function of the distance between islands. 
We have drawn two curves obtained by simulations for different heights of the wetting layer (z = 



and z = 10), and the curve obtained from Eq. IC without any adjustable parameter. The vertical 
dashed line corresponds to the two islands touching. Values of n^uf, used for the calculation of the 
curve are Usub = 240 for z = and Ugub = 90 for z = 10. 
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FIG. 4. Variation of Energy per particle as a function of the volume for a given value of the 
coverage 6. For each curve, only one point out of 50 has been precised in an aim to keep the figure 
clear. 
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FIG. 5. Variation of the volume of the islands and of the distance between islands at the 
minimum of energy as a function of the coverage 9. 
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FIG. 7. Variation of Energy per particle as a function of the volume for a given value of the 



coverage 6 with S{z) given by Eq. 18 and |T^. For each curve, only one point out of 50 has been 
precised in an aim to keep the figure clear. 
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